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ABSTRACT 

Simulations of Titan's atmospheric transmission and surface reflectivity have 
been developed in order to estimate how Titan's atmosphere and surface properties 
could affect performances of the Cassini radar experiment. In this paper we present a 
selection of models for Titan's haze, vertical rain distribution, and surface composition 
implemented in our simulations. We collected dielectric constant values for the Cassini 
radar wavelength (~ 2.2 cm) for materials of interest for Titan: liquid methane, liquid 
mixture of methane-ethane, water ice and light hydrocarbon ices. Due to the lack of 
permittivity values for Titan's haze particles in the microwave range, we performed di- 
electric constant (e r ) measurements around 2.2 cm on tholins synthesized in laboratory. 
We obtained a real part of e r in the range of 2-2.5 and a loss tangent between 1CT 3 
and 5.1CT 2 . By combining aerosol distribution models (with hypothetical condensa- 
tion at low altitudes) to surface models, we find the following results: (1) Aerosol-only 
atmospheres should cause no loss and are essentially transparent for Cassini radar, as 
expected by former analysis. (2) However, if clouds are present, some atmospheric 
models generate significant attenuation that can reach —50 dB, well below the sensi- 
tivity threshold of the receiver. In such 13.78 GHz radar would not be able 
to measure echoes coming from the surface. We thus warn about possible risks of 
misinterpretation if a "wet atmosphere" is not taken into account. (3) Rough surface 
scattering leads to a typical response of ~ — 17 dB. These results will have important 
implications on future Cassini radar data analysis. 

Key words: Titan; satellites, atmospheres; surfaces, satellite; radar 
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1 Introduction 



Titan is the second largest satellite in the solar system. It is the only one known to 
possess an extended and dense atmosphere, hiding its surface. Primarily composed 
of N 2 with less than 10% CH 4 and 1% H 2 , Titan's thick atmosphere is subject to an 
active chemistry induced by solar insolation, solar wind and charged particles from 
the saturnian magnetosphere. This chemistry results in the production of a great 
diversity of hydrocarb ons and nitriles in the stratosphere with increasing complexit y 



and molecular weight (jYung et al. 



1984! ; 



Toublanc et al. 



1995 



Lara et al. 



1996 



199 



Settling in the atmosphere, chemistry end products form an extensive system of organic 
aerosol haze responsible for Titan's characteristic orange color and hiding the lower 
atmosphere and the surface at UV-visible-NIR wavelengths. It is also possible that 
lighter hydrocarbons or nitriles diffuse downwards and condense as liquids or solids on 
aerosols below altitudes of 85 km. 

Ground-based observations, as well as Pioneer and Voyager missions, have brought 
a great number of constraints on the haze properties but revealed the difficulty to pre- 
cisely determine not only its horizontal and vertical structure but also its chemical 
composition. Concerning the surface, radar observations and microwave radiometry 
or near-infrared studies only reveal the presence of a heterogeneous surface plausibly 
compo sed of randomly distributed patches of water ice and organi c compounds, solid or 



liquid (IMuhleman et al. 



1990 



Coustenis et al. 



1995 



1997 



200 ll ). The precise nature 



of Titan's surface is still completely unknown. With the hope of disclosing the re- 
maining questions about its composition, the CASSINI-HUYGENS spacecraft - a joint 
NASA/ESA/ASI mission - was launched in 1997 and will reach Saturn in 2004. The 
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spacecraft, composed of an orbiter and a probe, will favour the exploration of Titan. A 
multipurpose radar was included in the scientific payload of the Cassini orbiter. Three 
active modes (SAR - for Synthetic Aperture Radar -, altimeter and scatterometer) and 
one passive mode (radiometer) will be used within the Radar experiment. Thanks to 
its active modes, the Radar instrument is expected to be able to observe Titan's surface 
through its dense atmosphere and thus determine whether oceans are present and map 
the geology and the topography of the solid surface of Titan. 

These scientific goals rely on analysis of the received radar echoes and the analysis 
depends on models describing microwave transmission in the atmosphere and on surface 
reflectivity. In this paper, we present a quantitative study of the possible impact of 
particles present in Titan's atmosphere on the Cassini Radar experiment and some 
results about surface reflectivity. We first present the models of Titan's surface and 
atmosphere we selected as a basis for our computations. Then, in section [3] we show 
results of experiments designed to determine the electromagnetic properties of titanian 
aerosol analogs - called tholins. We finally discuss in section H] the possible behaviour 
of the Cassini Radar in presence of atmospheric haze and rain on the basis of simple 
simulation of radar wave propagation for various models for aerosol and rain particles 
formation. We also show some simulations of the backscattered power by smooth and 
rough surfaces covered with tholins. 
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2 Distribution of condensates in Titan's atmosphere 
and on its surface 

2.1 Meso- and stratosphere: vertical distribution of aerosol 
haze 



Since 1973, there is observational evidence for the existence of an aerosol component 
in Titan's upper atmosphere. Particles with size greater than molecules were invoked 



to reconcile models o 



observations (jVeverka 



Titan's upper atmosphere and polarimetri c and photometric 



1973 



Zellner 



1973 



Rages and Pollack 



19801 ). Several radiative 



transfer models incorporating aerosols were published, considering the shape and op- 
tical properties of these aerosols as free parameters, as well as their vertical distribu- 
tion in terms of size and number density. These models were regularly adjusted to 



be consistent with the g rowing number o 



First generation models (IDanielson et al 



observational data of increasing r esolution 



1973 



Podolak and Danielson 



19771 ) assumed 



a homogeneous atmos phere where dust p a rticles and gas were uniformly m i xed above 



the surface of Titan. 



1980) , and 



Podolak and Giver! (1 19 79( 1. 



Rages and Pollack! ( 



Tomasko and Smith! ( 



Rages and Pollack! fll980h . 



Tomasko 



19831 ) used in their calculations multi-layer models and 
19821 ) used the first totally inhomogeneous models with a gradi- 
ent in particle size with the optical depth. They constrained the aerosol mean radius 
to between 0.1 fim and 0.5 fim, but always failed to simultaneously fit the polarimet- 
ric and photometric data acquired by Pioneer and Voyager missions when considering 
spherical particles of single size mode. 



Further studies greatly improved the model fit to data wit h the use o : 



complex haze distribution models. Inspired by the pioneer work of iToon et al. 



more 



(119801 ) 
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McKay et all (Il989l ) considered Titan's aerosols properties in the new light of micro- 
physical modeling, with the aim of integrating these results in a radiative-convective 
model of the thermal structure of Titan's atmosphere. They assumed that all particles 
at a fixed altitude had the same radius and that only thermal coagulation, sedimenta- 
tion, and electrostatic forces are responsible for variations in particle size and concen- 
tration with altitude. They predicted a size of ~ 0.4 /xm for the aerosols around an 
altitude of 85 km, in good agreement with previous values obtained from observational 



fits. In the same way, 



Cabane et al. 



(119921 ) developed a model including sedimentation 



and eddy diffusion accounting for transport processes, gravitational (coalescence) and 
thermal coagulation for collection processes. Their calculations focused on the pre- 
diction of the distribution and radius of aerosol particles throughout the atmosphere 
between the higher formation region and the lower stratosphere (above 85 km), be- 
fore condensation of hydrocarbons and nitriles on the aerosols could occur. The final 
particle rad i us ob tained in the low stratosphere is about 0.2 fim. In the meantime, 



Frere et al 



(Il990l ) pursued the development of a complete microphysical modeling of 



Titan's aerosols, including for the first time condensation processes of organics and 
nitriles prese nt in the gas phase i n the low str a tosph ere. Above 100 km, the pro files 



computed by 



McKay et al. 



(1989) 



Frere et al. 



fll990f ). and 



Cabane et al. 



fll992h are 



quite similar. The rece nt microphysical modeling taking into ac c ount aerosol growth 



into fractal aggregates (IRannou et al. 



1993 



Cabane et al. 



1993: 



Rannou et al. 

I r 



1995) 



gives concentrations and effective profiles (Fig. 1 and 3 in ICabane et al. I (119931 )) 



close to those given by classical microphysics dealing with spherical particles (close 
enough for the precision in radius and concentration we need), so we will not con- 
sider this model here. Besides in the wavelengths range we consider, fractal nature of 
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micron-sized scatterers would not have significant effect on scattering and absorption 
calculations. Figure CD- (a) shows profile s of aerosol r a dius a nd concentra t ion v ersus 



altitude resulting fro m computation by 



Cabane et al. 



fll992h . 



McKay et al. 



(jl989h . 



Frere et al. 



( 1199nh and 



2.2 Low stratosphere and troposphere: models for clouds and 
rain 



After a rapid lapse rate below 120 km altitude, the thermal profile of Titan reaches a 
minimum of 74 K at the tropopause (42 km). This thermodynamical cold trap region 
certainly induces saturation (even supersaturation for the dominant and more volatile 
species) and probably gas/liquid or solid phase transitions. These transitions should 
occur for most of the organics and nitriles pres ent in gas phase if seed nuclei are present 



( pagan and Thompson 



1984 



Frere et al. 



19901 ) . If nucleation can be activated, the most 



relevant scenario is that individual aerosol particles may act as condensation nuclei, 
accrete hydrocarbons and nitriles with more or less efficiency depending on solubility 

they pass through saturated regions of 



aeroso 


s, and 


( Lorenz 


1993a 



McKay et al. 



20011 ). However, observational evidence 



for a dense, persisten t and ubiquitous cloud co verage in the supersaturated regions 



is not yet conclusive. IGriffith et all (11998 



of Titan suggested a cloud c over. Titan's 



observed in September 1995 ( iGriffith et al. 



20001 ) hemispheric integrated observations 
arge cloud system that might have been 



19981 ) displayed cloud tops at an altitude 



of z ~ 15 ± 10 km, covering about 10% of Titan' s disk, and seemed to live for 



more than 2 days. Analysis of more recent spectra ([Griffith et al. 



2000) suggested 



another identification of tropospheric clouds. It revealed that clouds should reside at 
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the altitude of z ~ 27 km. These clouds covered ~ 0.5 ± 1% of the moon's disk, and 
dissipated in only 2 hours. Thus, if clouds exist, they seem to be patchy, occupying 
only a fraction of Titan's disc. It may also be an intense but transient and fickle 
phenomena, making them even more difficult to observe and firmly identify. 



Se veral models treating condensation in Titan's saturation 



veloped (IToon et al. 



1988 



Frere et al. 



1990 



Courtin et al 



9971) . aiming at the determination of cloud properties. 



1995 



region have been de- 



Samue 



Toon et al. 



dl988h and 



son and Mavo 



Courtin et al. 



(119951 ) studied the characteristics of hypothetic al pure met h ane c louds, considering at- 



Toon et al. 



(119881 ) assumed that if the 



mospheric radiative properties. In particular, 
mixing ratio of methane exceedzd 2% by volume near and below the tropopause, it 
would reach (and even exceed) saturation level, triggerin g its probab l e con densation. 



Toon et al. 



(119881 ) obtained 



Re-examining the Voyager IRIS infrared spectra of Titan, 
an optimal fit when adding a methane cloud composed of droplets with a radius be- 
tween 50 \im and 3 mm and a density between 1 and 3500 m~ 3 . This cloud cover, that 
can be considered as an ubiquitous homogeneous layer in addition to the aerosol haze, 



should extenc 



Courtin et al 



rom 10 km to 30 km. The new re-interpretation of IRIS spectra led 



(1995) to quite different conclusions. Their results, in best agreement 
with the data, did not contain a cloud layer. They found however that if a condensed 
CH4 cloud is added, it must have a mean particles radius of about 50 — 60 fim and 
be located near the tropopause (about 40 km in altitude), though they could not put 
constraints on its vertical extent. A combination of cloud characteristics inferred from 



Toon et al. 



(119881) and 



Courtin et al 



(119951 ) are summarized in Table 1. 



[Table 1] 
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Frere et al\ (119901 ) was the first to present a microphysical approach for conden- 
sation processes in Titan's low atmosphere, resulting in an estimation of the vertical 
radius and density distribution of aerosols. They elaborated a model dedicated to 
the determination of the distribution and chemical composition of the particles in the 
low atmosphere by taking into account differential condensation of several atmospheric 
compounds on aerosols in the cold trap region, in addition to classical transportation 
and collection processes. In their model condensation nuclei swell up, progressively 
covered with concentric layers of volatile compounds as soon as they reach saturation; 
supersaturation is not allowed. They constrained the size of particles by the volume of 
material available for condensation. They found that no significant condensation occurs 
above 100 km. The gas to liquid or solid phase transition starts modifying the particle 
structure below 90 km. From 0.1 to 1 //m around 100 km, the particle mean radius 
increases rapidly due to sustained accretion and reaches a few tens of microns at 85 km. 
At this stage of the condensation process, the major constituents surrounding aerosols 
are nitriles, mainly HCN and C^N 2 in solid state. A growing drizzle takes place down 
to the tropopause. Below 70 km, condensation of propane (C^Hg) becomes predomi- 
nant, then come C2 hydrocarbons (with decreasing altitude: C2H2, C<iH§ and C2H4). 
Between 30 km and the surface, the particle radius becomes greater than 200 fim, 
reaching a maximum of ~ 900 fim at 11 km. At these altitudes, liquid methane is 
the main constituent that condenses. Below 11 km, the strong positive temperature 
gradient triggers methane layer evaporation of the droplets. The fusion/sublimation 
rate is moderate until 2 km above the surface. Then, the process becomes very efficient 
and removes almost all the liquid around the aerosols. The particle radius decreases to 
about 50 — 100 fim. That gives rise to a mist of ethane-enriched aerosols just above the 
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surface (a "ghost" of the raincloud, consistent with considerations by lLorenzl (Il993al )). 
Variations of the droplet density follow those of the radius but in the opposite way due 
to colle ction processes. More recen t studies of rain occurence in Titan's l ower atmo- 



sphere ( jSamuelson and Mayo 



19971 ) led to similar conclusions as those of iFrere et al. 



(119901 ) . They obtained the same altitude for the maximum particle size (z ~ 11 



and ro ughly the sa me size pro fi 



( 119961 ) and recently 



Coll et al 



and siz e proposed by 



es ver sus altitude. The experimental work of 



12km) 



Frere et al 



McKav 



1999a) argued in favou r of drop l et ch emica 



ril990h (Fig. [B(b)). 



McKavi (|1996h and 



structure 



Coll et al. 



( I1999al ) demonstrated the insolubility of aerosol material in hydrocarbons and their 



relative solubility in nitriles. It is also con sistent with the identification of nitriles in 



solid phase in Voyager IRIS spect ra (CaNo (IKhanna et al. 



and HCN, HC 3 N, CH 3 CH 2 CN dCoustenis et al. 



1987 



Coustenis et al. 



19991) 



1999T ) that could constitute a 10//m- 



hail around the altitude of 90 km). Aerosol s might need a firs t nitrile layer to activate 



Frere et al. 



(Il990h (Fig. [B(b)). 



hydrocarbon condensation, as predicted by 
[Figure 1] 

2.3 Surface 

Very little is known about the surface of Titan, which is hidden from view by the organic 
haze. Voyager 1 and 2 flybys allowed some progress in the knowledge of the puzzling 
nature of the satellite's surface, though they failed to directly reveal surface details. 
Surface properties were only indirectly inferred from the atmosphere characterization 
performed by the spacecrafts. Voyager 1 and 2 flybys allowed the determination of the 
atmospheric properties with an accuracy such that it was possible to propose some new 
and more relevant models for Titan's surface nature and to postulate the existence of a 
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deep, global ocean of hydrocarbon (ILunine et a/.lll983l ). The fact that methane is likely 
to be depleted in the stratosphere within barely 10 million years due to a photolytical 
sink, coupled with the fact that Titan's surface temperature is close to the triple point 



of methane, supports the idea that Titan s 



rould possess a 



at its surface or at least in its subsurface ( Eshleman et al. 



iquid 



19831 ). 



hydrocarbon reservoir 



Observations have become available which suggest that Titan is not covered by 
an ubiquitous ocean. Subsequent ground-based radar observations, passive microwave 
radiometry and near-infrared observations give us today a new insight in the current un- 
derstandi ng of Titan's surface. D irect radar sensing of Titan's surface at a wavelength 



of 3.5 cm (IMuhleman et al. 



19901 ) revealed a surprisingly high reflectivity (much higher 
than the reflectivity of a global methane-ethane ocean) but still lower than that of the 
icy Galilean satellites. Radar reflectivity was also extremely variable from one night to 



the other, sug 
same time by 



gesting a heterogeneous surface. Passive radiometry undertaken at the 



Muhleman et al. 



( 119901 ) suggested an emissivity of about 0.9 at similar 
wavelengths (higher than those of Europa and Ganymede, but similar to that of Cal- 
listo). The only way they found to reconcile high emissivity and high reflectivity was to 
consider a surface covered with "dirty" ice (ice contaminated by organic compounds). 
This hypothesis was recently confirmed with the help of near-infrared observations from 
HST (Hubble Space Telescope), CFH (Canada- France-Hawaii) and, AEOS (Advanced 
Electro Optical System) telescopes in the methane-window regions (0.94, 1.08, 1.28, 
1.58, 2 and 5 jum), where it may be possible to probe the surface between methane 
absorption bands. These observations are consistent with d i rty ice, including areas 



of hydrocarbon seas (ILorenz 



1993b 



Lorenz and Lunine 



19971 ) 



Coustenis et al 



(2001) 



have investigated the nature of ice present at the surface of Titan. They took images 
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(albedo maps) of Titan's surface at 1.3 and 1.6 /im with a resolution high enough 



to resolve Titan's disk with 20 pixels. They obser ved the same 



equator t hat was previously o 



telescope ( IGibbard et al. 



Dserved by the HST (ISmith et al. 



brig 



it feature at the 



1996|) and the Keck I 



1999j). That allowed them to propose some candidates spec- 



trally relevant for the icy composition of the different (bright /dark) regions observed. 



By comparison to other surface spectra (ICoustenis et al. 



1995 : 



19971 ). they suggested 



the presence of methane (or ethane) frost, ammonia ice, and water ice on Titan's 
surface, probably c overed by sedimen ts of organics, as foreseen by photochemical at- 



mospheric models. iToon et al 



(119801 ) predicted indeed that the aerosols produced in 
the upper stratosphere should settle and deposit on the ground at the rate of 0.1 /im 
each year, to form a layer of about 400 m over geological times. The hydrocarbon and 
nitrile liquids they drag towards the surface are expected to stream down the relief and 
fill lakes or seas, whereas solid organic cores of the droplets should form sedimentary 
layers. It is nevertheless impossible to discriminate whether aerosols are intimately 
mixed with ice or whether they constitute a dark organic layer over an icy surface with 
some bright regions of exposed ice. The aerosols falling over a lake or a sea do not sink 
and create an organic layer, due to their poor solubility into hydrocarbon liquids. 



3 Electromagnetic characterization of titanian ma- 
terials 

3.1 Electromagnetic characterization of tholins 



Electric properties of titanian aerosol analogs were often measured in visible-UV range, 
but are poorly known in the microwave range. The only published value for dielectric 
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constant at GHz frequencies is the one taken from iKhare et al\ ( 119841 ) measurements. 
They found a permittivity of 4.71 — j0.013 at / ~ 300 GHz. Considering that tholins 
are major compounds of Titan's atmosphere susceptible to interact with Cassini radar 
electromagnetic wave, we performed a set of measurements of aerosols' analogs permit- 
tivity in the GHz frequency range. 

3.1.1 Synthesis of tholins 



Since the early work of IKhare and Saganl (119731 ). many laboratory simulations of Ti- 
tan's atmosphere have been achieved. Such experiments were conducted by irradiation 
with UV light, electrical discharge, or energetic electrons, of a gas mixture close to the 
Titan's atmosphere composition, isolated from contamination in a glass recipient - the 
reactor. Without any exceptions, these simulations produced a solid organic material 



- termed tholin - which h as optical properties sim i 



geometric albedo of Titan (IRages and Pollack 



1980 



ar to those needed to match the 



McKay et al. 



19891). 



In order to get tholin samples for further dielectric measurements, we conducted 
experiments in the LISA laboratory where a reactor essentially dedicated to the study 
of the titanian atmospheric chemistr y was designed. A det ailed descr i ption of the 



experiment al devices can 



De found in 



cording to 



McKay et al. 



Coll et al 



f|l999al ) and 



Coll et al 



(I1999bl ). Ac- 



(120011 ). this experiment promises to give the most accurate 
simulation of Titan's haze formation, as the dinitrile C^N? (detected in solid phase on 
Titan), has been identified for the first time as a product of t he simulation, this clearly 



constitutes a validation of the experiment (jColl et al. 



1999 



In our case, we performed three titanian atmosphere simulations with a low pres- 
sure flux (~ 1 —Ambar of total pressure) of a N 2 /CH 4: mixture (in a 98:2 proportion). 
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A glass tablet or a silica cylinder was inserted into the reactor where tholins are de- 
posited to form a film. In order to collect tholin samples, the reactor was installed 
inside a glove box filled with nitrogen to avoid contamination of the sample by the 
atmosphere in the room, in particular oxygen atoms, before electromagnetic charac- 
terization. Finally, the volume of deposited tholins - an essential parameter for future 
permittivity measurements - is determined using scanning electron microscopy to mea- 
sure the thickness of the sample. Table 2 summarizes synthesis conditions and volume 
estimations for each set of experiments. 

[Table 2] 

3.1.2 Dielectric constant measurements 

Tholin samples produced in the LISA reactor were characterized at the PIOM labo- 
ratory in order to measure their dielectric properties in the GHz range. Permittivity 
measurements were performed using a rectangular resonant cavity pierced on its top 
with a hole used to insert the sample. For the purpose of these experiments, we used 
two different kind of cavities : one at 10 GHz (i.e. A = 3 cm) and a second one at 2.45 
GHz (i.e. A = 12.2 cm). As the change in the dielectric constant value should be small 
in the range 1-10 GHz, the choice of a cavity was made only because of the sample 
holder shape : the 10 GHz cavity is too small to allow a large rectangular hole on its 
top center and only a cylindrical sample holder can be used. Each cavity is composed 
of a wave-guide short-circuited at each end. Excitation and detection of resonance 
inside the cavity are done by co-axial transitions mounted on connectors linked to a 
network analyzer (ANRITSU 37325A). 

The resonance modes that can be established in a cavity are deduced from the 
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theory of wave propagation along a wave-guide. These modes are characterized by 
three integer indices m, n and p. In the case of a rectangular cavity, TEOlp resonance 
modes are usually used (m = and n = 1 as dominant mode, p representing the 
number of maxima of the electric field along the wave-guide length). In particular, the 
index p is chosen odd (p = 3) in order to ensure the position of a maximum in the 
center of the cavity, where we insert the sample. The sample holder symmetry axes 
are set in a way that they are parallel to the electric field established in the cavity, in 
order to minimize perturbations due to depolarization effects. 

Once the cavity is connected to the network analyzer, we first introduce an empty 
sample holder into the cavity to measure the resonance frequency F , and the quality 
coefficient Qo of the cavity as blank references. Then, we insert an identical sample 
holder - same glass/silica volume - but covered with a tholin film. The introduction 
of a tholin sample triggers a deviation from the resonance frequency, AF, and an 
additional attenuation conducting to a new quality factor Q\. As it is well established 
for resonant cavity experiments, the measurement of AF, Q , and Qi allows to directly 



obtai n the dielectric constant e — je of the sample thanks to the relations ( Boudouris 



19711 ): 



AF . i 

— = -a(e - 1) (1) 



A(^- - = 2ae" (2) 

with a = «o x Vgampie, V sa ra V ie being the volume of the tholin sample and «o being 
the cavity filling coefficient (ratio between the tholin sample + sample holder volume 
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and the cavity volume). 

For each set of tholin samples (see Table 2), we obtained a value for the tholin 
dielectric constant. Results of our measurements are shown in Table 3. While real 
part of the dielectric constant shows a very little variation (~ 10%), the imaginary 
part varies by a factor of 50. This variation is due to the poor estimation of factor 
Q, which is highly sensitive to the volume of material introduced in the cavity. An 
accurate value for this volume is very hard to obtain. 

[Table 3] 

3.2 Electromagnetic characterization of other titanian mate- 
rials in liquid and solid states 

Apart from tholin characterization, we also searched in the literature for values of 
permittivities of materials of interest for radar sounding of Titan's atmosphere and 
surface, in the microwave range. 

3.2.1 Pure liquid methane 

The dielectric proper ties of pure liquid methane and their implications for Titan's study 



were investigated by 



Thompson and Squyresl (Il990l ). They published a rough estimate 



of e r for liquid CH4 around 94 K: s ~ 1.7 and e" ~ 1.5 x 10~ 2 . This is fully consistent 
with the value found in Handbook of Chemistry and Physics (1988-1989) around 100 K 
in the hundreds of MHz range. 
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3.2.2 Titan's raindrops and hypothetical surface oceans: liquid hydrocar- 
bon and nitrile mixtures 



According to iFrere et all ( 119901 ) . Titan's raindrops are primarily made of light hy- 
drocarbon liquids instead of pure methane, progressively accreted in concentric layers 
around aerosols. During the fall of the raindrops, internal dynamics stirs the liquid 
CH^ and C%H.§ into a homogeneous mixture. Such a mixture was th e subject of a few 



experi mental studies aiming at determining its dielectric properties. 



Muhleman et al. 



(119901 ) indicate that the real part of the dielectric constant for liquid hydrocarbons is 



well known to be in the range of 1.6 to 2 and t hat lo ss tangents are of the order of 



10 3 — 10 4 in the GHz range. 



Singh and Miller! (119791 ) found a real part of permittiv- 



ity between 1.77 a nd 1.9 for methane-e thane mixtures at low temperature (typically 



around 100 K) and 



Ulamec et al 



(1992) o btaine d values between 1.71 and 1.95 around 



94 K. The values from ISingh and Miller! (11 9791 ) were obtained using a fit to an equa- 



tion derived from the Clausius-Mossott i function, on t he ba sis of liquid compressibility 



considerations, and those coming from 



Ulamec et al. 



(119921 ) were inferred from capac- 



ity measurements. We have no indications about what could be frequency range they 



considered. E 
of 



Sen et al 



owever, such values are consistent with dielectric constant measurements 



19921 ) for liquid alkanes and hydrocarbon mixtures at 1.2 GHz and 298 



K. 



3.2.3 Titan's surface: hydrocarbons and water ice 



Titan's surface is supposed to be composed of various icy materials (H?Q ice, hydro 



carbons in solid state) due to its low 94 K temperature (ILellouch et al. 



19901 ) . probably 



fractionally covered or mixed with tholin material. It then should present a high and 
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contrasted reflectivity in the radar wavelength range. Dielectric constants of such ma- 
terials were experimentally measured. We had access to precise values for the real part 
of their permittivity, but only to rough estimations for the imaginary part: H2O ice cor- 
responds to £ ~ 3 .1, measured aroun d 10 GHz and valid within the temperature range 



from -1 to -185°C f Ulabv et al. 



(119821 ) . vol. 3, appendix E, p. 2 028), and a C^-C^b 



ice m ixture has e ~ 2 — 2.4 in the same wavelength range ([Thompson and Squyres 



1990). 



All the dielectric constants of interest for radar sounding of Titan are summarized 
in Table 4. 

[Table 4] 

4 Simulation of the CASSINI radar behaviour 
4.1 Testing the radar penetration through the atmosphere 

In this part, we present simulations of the CASSINI Radar experiment in active modes. 
This section mainly deals with the dielectric models we have developed to quantify the 
impact of condensates in Titan's atmosphere. 

We used in our calculations the complex dielectric constant we measured for 
tholin material. Uncertainties on the permittivity of tholins have repercussions on 
the attenuation we calculate: concerning atmospheric transmission simulations, the 
attenuation does not vary linearly with the dielectric constant. The minimum value 
for attenuation was obtained with s =2.1 and e = 0, and the maximum attenuation 
corresponds to e — 2.1 and e = 0.1. For the preliminary results on homogeneous 
atmosphere models, we considered the mean value measured for the dielectric constant 
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of tholins: e r = 2.2 -j'0.05. 

4.1.1 Description of the Rayleigh-Mie scattering model 



We consider the propagation of a radar wave through Titan's atmosphere. The simple 
model we developed enables us to determine the attenuation an incident radar wave 
(emitted power p mc ) will undergo while crossing the dense atmosphere of Titan. The 
attenuation in decibels is given by 

pground 

Att tota i(dB) = Att gases + Attpartides = 10 log w { p inc ) (3) 
with paround bging the rac l ar wa ve power reaching the surface of Titan. 

The gaseous component is considered with no ambiguity as totally transparent 
for the radar and its permittivity is set equal to one. Thus, 



Att gases (dB) = (4) 

Only the condensate - aerosols, hail and rain particles - component could have 
sufficient size, density and permittivity to interfere with the radar wave. In this case, we 
can reasonably simulate the medium that the wave will encounter before reaching the 
surface by a diluted medium composed of dielectric spheres in suspension in vacuum. 
Only single scattering and absorption of the incident wave are considered here and we 
can then write for an homogeneous atmosphere: 



pground 

— e -[C(o- S cat+a a bs)-^e)] (5} 

pine \ ' 
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with, 



C, particle density, 

h, thickness of the atmosphere, 

9, incidence angle of radar wave (set to 0° for the altimetric mode), 
a sca t, single scattering cross section, 
<^abs, absorption cross section. 

Cross sections for spherical scatterers are taken from Rayleigh and Mie theories. 
Transition between the Rayleigh and Mie regimes occurs when the particle radius 
is greater than r^ m = 0.05. Rayleigh formalism is valid only for small particles in 
comparison to the wavelength. In the case of the Cassini Radar, the upper limit to the 
radius ru m of the scatterer is set to 1.08 mm. 

Raylei gh scatter i ng an d absorption cross sections can be analytically solved and 



are given by 



Ishimarul (119971 ): 



1287r 5 r 6 ,£ 



" scat 



3A 4 1 e r + 2 1 



(6) 



lir 2 r 3 e , 3 



Cabs 



3A 



(7) 



Mie scattering cross se ctions are numerically comp uted by means of a clas- 
sic BHMIE code taken from iBohren and Huffman! (119831 ). The complex parameter 



e r = s — je" represents the dielectric constant of the spheres of radius r. It should 
be frequency-dependent and was experimentally measured around the Cassini Radar 
experiment frequency (cf. section ED- 

[Figure 2] 

In our model (for heterogeneous atmosphere cases), the atmosphere of Titan is 
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vertically cut off into homogeneous layers containing particles of fixed size and density 
(Fig. [2]). We then compute the ratio between outgoing and incident flux P° ut / P\ nc for 
each layer % of thickness hi. The resulting final attenuation is the sum of attenuations 
induced by each layer i. It is given by : 



Equation (jSJ) indicates that the part of the flux scattered or absorbed by the 
particles is lost. It reproduces well what could happen to a wave in single scattering 
regime along the radar line-of-sight. 

4.1.2 Homogeneous atmospheres 

We computed the attenuation of the centimetric signal emitted by the Cassini Radar 
instrument induced by the presence of particles in a ideal Titan-like homogeneous 
atmosphere. In this simple case, we considered single crossing of the wave through a 
totally homogeneous atmosphere, made of a unique layer with a single mode of particle 
size, density and dielectric properties. 

The thickness H of the homogeneous layer was constrained by the aerosol haze 
extension within Titan's atmosphere and was fixed to h = 500 km. We fixed e r to its 
experimental mean value of 2.2 — j0.05. The only free parameters in our simulation 
were particle radius and density. We decided to put a constraint on these parameters 
from extreme values found in literature: density ranges between unity and a maximum 
around 10000 cm -3 , and radius ranges between a few nanometers and a high value of 
20 fim. Attenuations induced by such particles over a 500 km thick layer are shown 
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in Fig. (a). Each point of the figure represents the attenuation calculated for a 
unique atmospheric scenario, ie a unique pair of particle radius and density over a 
homogeneous Titan atmosphere. 

[Figure 3] 

The detection threshol d (noise equivalent o"o) of the Cassini Radar instrument 



is limited to about —25 dB (jElachi et al. 



199ll ). It means that, with the hypothesis 



that all the power reaching the surface of Titan is backscattered to the radar, a signal 
attenuation greater than 12.5 dB for the first crossing of the atmosphere would lead 
to no signal detection by the instrument. Figure [3j-(b) is a 2-D representation of the 
variation of the attenuation versus particle radius and density. It shows the position 
of the 12.5 dB attenuation limit, and attenuation greater than 12.5 dB for only one 
atmosphere crossing is represented in gray on the Fig. [3J-(b)). Some radius/density 
combinations are then prejudicial to the Radar experiment and its main objective 
of surface imaging. The model we considered here is very simple but shows that 
it is worthwhile to study the possible impact of atmospheric particles (aerosols, hail 
particles, drizzle and rain droplets) on the radar experiment. It allows also to define 
some first order ranges of radius and density for which the radar could encounter 
difficulties to see the surface of Titan. It would be the case with, for example, a 
500 km thick atmosphere constituted of particles of 6 /im radius with a 1000 cm -3 
density, which is not completely unexpected or unrealistic at first order when we refer 
to distributions retrieved from observations or microphysical calculations including rain 
occurrence ( particles are then concentrated in a thinner layer, but are much bigger than 



ten microns (IToon et al. 



1988 



Frere et al. 



1990l )). The scattering cross section increases 



23 



with the sixth power of the radius of the scatterers and constitutes the more sensitive 
parameter. These first results lead us to the conclusion that Titan's atmosphere might 
not be as transparent to Cassini radar frequencies as it was expected. Radar pulses 
could at least be affected and distorted, even if they encounter smaller particles than 
those we considered in this section. 

4.1.3 Heterogeneous atmospheres and cloud layers 



In this section, we take into account as a modeling refinement the possible heterogene- 
ity of aerosol properties with altitude. We do not treat the atmosphere as an ideal 
homogeneous layer anymore but as a more complex medium as predicted by numerous 
published models. 

We used the heterogeneous distribution of particles discussed in section |2] in order 
to simulate a more realistic Titan's atmosphere (see Fig{T]-(a)). We also included as a 



reference two semi- homogen eous models for T i tan's atmosphere inferred 



and Voyager image analysis (jSmith et al. 



1981 



1982 : 



Tomasko and Smith 



rom Pioneer 



19821 ) . High- 



resolution views of the limb at high-phase angles revealed the presence of several layers 
of aerosols: an optically thin layer of haze ~ 50 km thick roughly 100 km above a main 
aerosol layer. The aerosol particle average size in the main haze layer was found to be 
between 0.1 and 0.5 \im for a real part of refractive index n r chosen between 1.7 and 2 



ue for the partic 



Toon et al 



e radius 



fll980h . They 



at visible wavelengths. For the upper haze layer, a mean val 
was taken to 0.3 fim. Particle densities were calculated by 
found 10 4 cm -3 for the first tenth of the visible optical depth in Titan's atmosphere 
and 10 3 c m" 3 at an optical d epth of 1. In addition to the detached haze layer, and 



above it, 



Rages and Pollack! (119831 ) discovered a small but persistent feature in the 
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Voyager 2 extinction profiles, characteristic of another haze layer with the following 
physical properties: roughly 50 km thick, composed of 3 /wn-particles with a density of 
~ 0.2 cm~ 3 (cf. Table 5). We considered a mean aerosol permittivity e r = 2.2 — j0.05. 
Errorbars on attenuation come from uncertainties on tholin permittivity (the minimum 
value for attenuation was calculated with e =2.1 and e" = 0, and the value maximum 
with e = 2.1 and e" = 0.1). 

We considered also the probable occurrence of hydrocarbon rain at low altitudes. 
It was treated differently according to two sce narios: (1) We sim ula te the rain as 



anoth er independent layer in addition to aerosols (IToon et al. 



fll988h and 



Courtin et al. 



(119951 ) a pproach) an c 



aerosols (IFrere et al 



(2) we consider rain as the result of continuous growing of 



(119901 ) approach). In case (1), we construct the cloud cover by 



substituting the "dry" particle distribution (two homogeneous multi-layer models as 



defined earlier and models pr esented on 



methane droplets inspired by 



Toon et al. 



ig. [T -(a)) by a homogeneous layer containing 



(119881 ) (cf. Table 1). The vertical extension of 



the cloud is supposed to range between altitudes of 30 and 10 km. We fixed an average 
methane droplet radius to 2 mm and their permittivity to e r = 1.7 — j0.015, based on 
the hypothesis that the raindrops are only composed of pure liquid methane and that 
the tholin nuclei can be neglected (the liquid phase hugely dominates the composition of 
th e droplet). The way w e contr ol the density of the clo ud was inspired by recent work 



of 



Samuelson and Mayol (119971 ) and 



Guez et al. 



(119971 ) . They suggested that liquid 



hydrocarbons may condense only from a limited subset of the aerosols fa 



Guez et al 



ling through 



(119971 ) introduced 



the saturated layer, in particular for kinetic reasons, 
the concept of nucleation rate. Indeed, the condensation process efficiency is highly 
dependent on the contact angle of condensed materials on the nucleus. This is an 
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essential parameter to estimate the compatibility between nucleus and the condensed 
phase. As contact angles are very unconstrained parameters, we treated the problem 
in a simple way thanks to a condensation efficiency factor / that can be free to vary 
between and 100%. Using this coefficient, we simulate the fact that only a fraction of 
aerosols can supply the cloud in droplets at the condensation level: for an atmosphere 
sampled every kilometer, cloud concentration is then given by Cdro P iets( z = 30 km) = 
f x Caerosois(z = 31 km). We checked that densities obtai ned for the t o p of t he methane 



Toon et al. 



(11988( 1. whatever 



cloud (z = 30km) never exceed the limits imposed by 
higher aerosol distribution can be. The f parameter was fixed equal to 10 _4 %, giving 
particle concentrations at the cloud top of 1000, 1000, 8, and 250 m~ 3 respectively for 
the atmospheric models A, B, C and D (cf. Table 5). The cloud's particle radius and 
concentration are then constant with respect to the altitude from the cloud top (30 km) 
to its base (10 km). Aerosols which don't participate in the condensation process are 
neglected, and Titan's atmosphere is then only constituted of a homo g eneous cloud of 



Lorenzl (ll993al ) showed 



methane rain droplets between 30 and 10 km. Below 10 km, as 
pure methane droplets should evaporate and resupply aerosol layers, forming a "ghost" 
mist of aerosols close to the surface. In order to simulate rain sublimation in the last 
kilometers of Titan's troposphere, we stop the rain layer below 10 km and replace it 
by the "dry" scenario we used to describe the atmosphere above 30 km. 



In case (2), we simply used the heterogeneous distribution computed by lFrere et al. 



( 119901 ) (Fig. [Q- (a)) as the complete description of the fate of particles over the full ex- 
tent of Titan's atmosphere, even at low altitudes where condensation should occur. 
As they treated the differential condensation of a great number of hydrocarbons and 
nitriles, we chose an averaged permittivity for the raindrops inferred from those mea- 
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Frere et al. 



(119901 ) thanks to modeling of 



sured for a homogeneously mixed hydrocarbon-nitrile liquid: e r = 1.8 — j0.002 (as 
for case (1), electrical properties of the condensation nuclei were neglected). In this 
case, the artificial factor / has no sense anymore. All aerosols participate in the con- 
densation process. The density and the si ze of the ra i n dro plets below 100 km are 
thus self-determined by the calculations of 
collisional and sticking processes. 

All the models used in our calculations are summarized in Table 5. 

[Table 5] 

[Figure 4] 

Figure H] shows the results obtained for the nine models described in Table 5. 
Cases A, B, C and D correspond to attenuations induced by a dry atmosphere, only 
constituted of aerosol hazes. The attenuation never exceeds -0.02 dB, even for the least 
favorable scenario B (containing the highest particle size and density over the greatest 
part of the atmosphere), and are globally negligible (< — 10~ 3 dB). We can conclude 
that, without rain occurrence in the troposphere, Titan's atmosphere should be totally 
transparent for the 13.78 GHz radar of the Cassini orb iter. When we consider cases 
Ajj, Bjj, Cjj, DJj and E, the conclusions are quite different. The addition of a rain 
layer in the middle troposphere increases the signal attenuation by at least four orders 
of magnitude. For DJj, Ajj-BJj (rain attenuation overflows the attenuation induced by 
aerosols and the difference between scenario Ajj and Bjj is then small; thus Ajj and Bjj can 
be considered as one unique scenario), and E models, attenuation reaches respectively 
-5 dB, -20 dB and -25 dB. On the other hand, the value of attenuation is still low for 
scenario Cjj (-0.17 dB), due to the fact that the rain cloud modeled for this case presents 
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a very low density. This parameter is controlled by the factor /, fixed to 10 _4 %, and 
the density of aerosols just above 30 km. For scenario CJj, the cloud has a density 
of only 8 m -3 , in comparison to models Ajj and DJj which have a raincloud density 
of respectively 1000 m~ 3 and 250 m~ 3 . We have to keep in mind that the titanian 
aerosol modeling of 



McKay et al. 



(119891 ) (scenario C and CJj) was rather simple, and 



the aerosol de nsity it predicts a round 30 km is ce rtainly less accurate than the one 



obtained from 



Frere et al 



Cabane et al. 



( 119921 ) (scenario D 



(119901 ) (scenario E) and 
and DJJ). We are convinced that, even if clouds in Titan's atmosphere are transient and 
rare phenomena, they could surely affect the radar signal. For two of the scenario we 
studied (Ajj and E), clouds should totally attenuate the radar pulse before it reaches 
Titan's surface. For model DJj, it should also cause sufficient attenuation to possibly 
blur surface images by diminishing the contrast between distinct surface echoes. 

[Figure 5] 

The rain layer we used in addition to models A, B, C, and D in order to simu- 
late the condensation process in Titan's low atmosphere is intrinsically artificial in its 
construction. We explored the way the attenuation may vary if we free some badly 
constrained parameters of this layer: radius of droplets and condensation efficiency f. 
The rad ius can vary b etween mm (no rain) and 4.5 mm (limit imposed by aerody- 



namics (jLoren2 



1993al )). For the parameter /, even if there are no means at the present 



time to physically impose a constraint on it, we can consider that the value we chose 
(/ = 10~ 4 %, i.e. one aerosol per millio n acts as a conde nsation nucleus) in order to 



match the cloud densities predicted by 



Toon et al. 



(119881 ) is very low. Higher values 



(up to 0.1%) cannot be totally excluded and should cause much more attenuation. We 
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present in Fig. [5] the results of this parametric study: the attenuation as a function of 
rain droplet radius (Fig. O-(a)) and condensation efficiency (Fig. O-(b)) for models Ajj- 
BJj, Cfj and DJj. The increase of r and / induces a drastic rise of attenuation for all the 
cases considered. Attenuations are much greater than unity as soon as r reaches 2 mm, 
and / reaches 0.001%. Figure [5] demonstrates the high sensitivity of the attenuation 
to r and /, which are unfortunately weakly constrained. 

4.2 Surface reflectivity 

As aerosols are supposed to deposit on Titan's surface and form a "dirty ice" mixture, 
we tried to simulate the radar response of a surface covered with tholins. It allows us 
to refine the first hypothesis we made in the previous section of a perfectly reflecting 
surface (i.e. all the energy reaching Titan's surface is backscattered to the instrument). 



Since the Cassini Radar experiment wil 



A = 2.17 cm, HH polarization (jElachi et al 



to appear as rather rough for the radar ( 



oper ate in the Ku-band (/ = 13.8 GH z, 



1991), most of Titan's surface is likely 



Dierking 



1999). We shall then consider 



surface scattering models developed for medium rough to very rough surfaces, i.e. 
the Physical Optics (PO) model and the Geometric Optics (GO) model, which were 
derived from the Kirchhoff model under the scalar approximation and the stationary- 
phase approximation respectively. These models are relevant to evaluate performances 
of a side-looking imaging radar. We shall consider here a single homogeneous layer 
covered by tholins (characterized by its averaged dielectric constant e = 2.2 — 0.05j), 
water ice or liquid methane, and characterized by its surface roughness. The latter is 
defined by the surface correlation function p(x) (a Gaussian function is assumed here), 
the height standard deviation a, and the correlation length L (cf. Fig. E]). 
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[Figure 6] 

The PO model validity range is defined by the following conditions: 



V2a 



< 0.25, kL > 6 and L > 2.76<xA 



(9) 



where k is the wave number of the incident plane wave (k = 2tt/\). The HH non- 
coherent sca t tering coefficient is then given as a function of the incidence angle 9 by 



Ulaby et al. 



(11982h : 



a° HH (9) = 2k 2 cos 2 6Y H {e)e-^ k,y cos6 ^ ■ V(4fc 2 cr 2 cos 2 6) n /nl / p n (x) J (2kxsm9)xdx (10) 

, Jo 



where Jo is the zeroth-order Bessel function of the first kind, Th(0) —\ Rh{6) | 2 is 
the Fresnel reflectivity, and p(x) = exp(— x 2 / L 2 ) is the Gaussian surface correlation 
function. 

The GO model validity range is expressed by: 



(2kacos6) 2 > 10, kL > 6 and L 2 > 2.76a X 



'ID 



and the HH (or VV) noncoherent sca t tering coefficient is given as a function of the 



angle of incidence 9 by 



Fung and Eoml (Il98ll ) 



r(o)e- 



-tan 2 (6»)/2m 2 



12) 



2m 2 cos 4 6 

where m = is the rms slope for a Gaussian surface and T(0) is the Fresnel reflec- 
tivity evaluated at normal incidence: 



30 



1 — Je/u,- 

r(o) = | — ^==| 2 (13) 



1 + yj£/fi 

(we consider here a magnetic permeability \i — 1). 
[Figure 7] 

Figure [7] presents the radar backscattered power for two types of surfaces: a 
relatively smooth surface with a = 0.5 cm and L = 25 cm described by the PO model, 
and a rough surface with a = 7 cm and L — 10 cm described by the GO model. 
Calculations were performed for surfaces covered by tholins, water ices and liquid 
methane. Uncertainties on experimental tholin dielectric constant were considered, 
but the resulting imprecision on backscattered power actually stays between the values 
obtained for water ice and liquid methane. For the mean look an gle foreseen for the 



199ll )). we can see that 



imaging mode of the Cassini radar (8i ~ 11° ( lElachi et al. 
the backscattered power for a smooth flat surface is very low (about -80 dB), much 
l ower than the noi se equivalent ctq of the instrument announced to be around -25 dB 



(lElachi et al. 



19911 ) . On the contrary, a rough and flat surface should present a much 
higher return, about -17 dB. As a first approximation, the Cassini imaging radar should 
only see very rough regions covered by tholins (land surfaces ?) whereas smooth and 
flat surfaces (seas and lakes ?) should appear as dark zones. Nevertheless, for non-flat 
smooth regions, Fig. [7] shows that slopes facing the radar with an incidence angle 9 
less than 7° would produce a backscattered power higher than -25 dB. Considering 
the mean look angle 8l, it corresponds to a slope ranging between 4° and 28°, which 
represent a large fraction of slope values that can be observed on planetary surfaces. 
Results obtained here should be of course considered as indicative, since very little is 
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known about Titan's surface composition and roughness. In particular, due to the low 
loss tangent of tholins (tan 5 = e je = 0.02) and other surface materials such as water 
ice, the penetration depth Dp of the incident Ku-band wave given by: 

^ = S { fV 1 + ( 7 )2 - irI/2 (14) 

is about 10 cm. This could give rise to some attenuation and volume scattering 
effects that could change the radar backscattered power computed from Eq. ffTUj) and 
( I12p . A mixture between silicate and tholins should also lead to higher values of e, and 
then to a higher reflectivity of the surface. 

If we make the link between the results we obtained for atmospheric attenuation 
and the surface backscattered power, it leads to more realistic conclusions for the 
attenuating role of the atmosphere. The surface should not be perfectly reflective 
for the radar wave as we hypothesized in the previous section. In the favorable case 
of a rough surface, the reflected signal could still lose 17 dB, and the limit of the 
acceptable atmospheric attenuation should then be reduced to -4 dB instead of -12.5 
dB if considering a radar sensitivity limit of —25 dB. This dramatically increases the 
impact of an attenuating atmosphere on the Radar experiment (see Fig. @] and [5] with 
the attenuation limit fixed at -4 dB). 

5 Conclusion and perspectives 

The Cassini mission has been designed to send a spacecraft to the planet Saturn, and 
deploy an instrument probe, Huygens, that will descend to the surface of Saturn's moon 
Titan. The radar instrument onboard the orbiter is expected to bring the scientific 
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community new information to characterize the surface of Titan. We are convinced that 
correct interpretation of the measurements made with the Cassini Radar instrument 
requires simulations and laboratory measurements to anticipate possible atmospheric 
effects and then better invert and interpret radar data. 

In this paper we presented a study of the behaviour of a 13.78GHz plane wave go- 
ing through Titan's atmosphere and backscattered by its surface. Titan's atmosphere 
was modeled by a diluted medium filled with spherical particles of size, density and 
permittivity estimated from what could be found in the literature and laboratory ex- 
periments. We considered one homogeneous atmosphere model and five heterogeneous 
models including or not condensation that could occur in the lowest part of titanian 
atmosphere. In order to fix some parameters of our models, we conducted a series of 
experiments in order to derive the dielectric properties around 13.78 GHz of Titan's 
aerosol analogs produced in laboratory. 

The results we obtained in terms of radar power attenuation after crossing the 
simulated Titan's atmosphere are of two kinds: 

• the "dry" heterogeneous atmospheric models, only taking into account the aerosol 
component, do not cause any attenuation and the radar wave will reach the 
surface of Titan without any losses. The particle radius never exceeds 1 \im (far 
from the critical limit around 10 \im we estimated from the ideal homogeneous 
case) and the corresponding attenuation never exceeds —0.01 dB. In such cases, 
Titan's atmosphere is totally transparent to microwaves and its effect should be 
totally neglected in the future processing of Cassini radar data. 

• on the contrary, when we add in our simulation a rain layer in the last few 
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kilometers, the attenuation reaches and rapidly exceeds the instrument sensitivity 
limit, due to the enhancement in size of the particles encountered (that could 
reach a few millimeters). A wave emerging from such a cloud layer could be so 
attenuated that the radar antenna wouldn't be able to detect a returned echo 
from Titan's surface. This could be a problem for the retrieval of Titan's surface 
images, as well as for altimetry measurements, and it should be taken into account 
when interpreting future data of the Cassini Radar experiment. 

Simulation of the backscattered signal by Titan's surface covered by tholins also shows 
that a rather low return could be expected, even for rough surfaces. This could dra- 
matically increase the possible attenuation role of Titan's atmosphere. 

There is also a non-negligible risk for the Cassini radar instrument to detect 
a false signal coming from a reflective atmospheric layer, screening and flooding the 
surface signal, that could cause false interpretation of "surface" images. As the Cassini 
radar is a multipurpose instrument, a way to discriminate between possible atmospheric 
and surface echoes would be to combine altimetric and SAR data. It could then be 
possible to avoid false interpretation of "non-surfacic" signal, if both altimetric and 
SAR acquisition can be performed in the same short time period. The altimetry 
pulse shape could also be analyzed in order to detect an atmospheric effect and get 
unexpected information about Titan's lower atmosphere (rain occurrence and extension 
of cloud systems, size, density and velocity of the particles). In a flyby strategy which 
would maximize the coverage of Titan's surface, the radar will certainly not map the 
same region twice. In this case, using the altimetry data would be the only way 
to discriminate possible cloud layers from surface features. On the other hand, two 
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flybys over the same area could give the opportunity to distinguish persistent features 
related to the surface from transient ones, which are likely to be due to atmospheric 
phenomena. 

The new analysis scheme we propose for the Cassini radar measurements will 
be subject to further investigations. For this purpose, we shall improve our present 
simulation of the Cassini Radar instrument in order to model the behaviour of the 
pulse it will send towards Titan's surface (not only power attenuation, but also shape 
distortion). We shall also take into account the doppler shift in frequency of the 
signal sent and received by the orbiter, as could be produced by the relative motion 
of the atmospheric scatterers. In order to properly prepare for the analysis of future 
Cassini radar data, we also need for more accurate measurements of haze and rain 
properties (size, density and dielectric constant) as we showed the sensitivity of our 
calculations to these parameters. Rain droplet radius and condensation efficiency in 
the low atmosphere are in particular crucial parameters. Such information could also 
be retrieved from the combined interpretation of independent observations of Titan's 
atmosphere coming from the other instruments onboard the Cassini orbiter or Huygens 
probe, as well as results of future rain modeling. 

ACKNOWLEDGMENTS 

We sincerely thank PIOM's director, Jean-Paul Parneix, and LISA's director, 
Frangois Raulin, who welcomed us in their laboratories and gave us the opportunity to 
instructively collaborate with some of their team members. We are also very grateful 
to Jonathan Braine who patiently checked the English language of the paper. This 
work was supported by the French Programme National de Planetologie of the Institut 



35 



National des Sciences de l'Univers, CNRS. 

References 

Bohren, C.F., and D.R. Huffman 1983. Absorption and Scattering of Light by Small 
Particles, A Wiley-Interscience Publication. 

Boudouris, G. 1971. Cavite electromagnetique, , Dunod, Paris edition. 

Cabane, M., E. Chassefiere, and G. Israel 1992. Formation and growth of photochemical 
aerosols in Titan's atmosphere. Icarus 96, 176-189. 

Cabane, M., P. Rannou, E. Chassefiere, and G. Israel 1993. Fractal aggregates in 
Titan's atmosphere. Planet. Space Sci. 41, 176-189. 

Coll, P., D. Coscia, N. Smith, M.-C. Gazeau, S.I. Ramirez, G. Cernogora, G. Israel, 
and F. Raulin 1999a. Experimental laboratory simulation of Titan's atmosphere: 
aerosols and gas phase. Planet. Space Sci. 47, 1331-1340. 

Coll, P., J.C. Guillemin, M.-C. Gazeau, and F. Raulin 1999b. Report and implications 
of the first observation of C 4 iV 2 in laboratory simulations of Titan's atmosphere. 
Planet. Space Sci. 47, 1433-1440. 

Courtin, R., D. Gautier, and CP. McKay 1995. Titan's thermal emission spectrum: 
reanalysis of the Voyager infrared measurements. Icarus 114, 144-162. 

Coustenis, A., E. Lellouch, J. P. Maillard, and CP. McKay 1995. Titan surface: Com- 
position and variability from its near- infrared albedo. Icarus 118, 87-104. 



36 



Coustenis, A., E. Lellouch, M. Combes, R. Wittemberg, CP. McKay, and J. P. Mail- 
lard 1997. Titan's atmosphere and surface from infrared spectroscopy and imaging. 
In Astronomical and Biochemical Origins and the Search for Life in the Universe 
(C.B. Cosmovici, S. Bowyer, and D. Werthimer, Eds.), 277-234. Proceedings of IAU 
Colloquium 161. 

Coustenis, A., B. Schmitt, R.K. Khanna, and F. Trotta 1999. Plausible condensates 
in Titan's stratosphere from Voyager infrared spectra. Planet. Space Sci. 47, 1305- 
1329. 

Coustenis, A., E. Gendron, O. Lai, J. P. Veran, J. Woillez, M. Combes, L. Vapillon, 
T. Fusco, L. Mugnier, and P. Rannou 2001. Images of Titan at 1.3 and 1.6 /im with 
adaptative optics at the CFHT. Icarus 154, 501-515. 

Danielson, R.E., J.J. Caldwell, and D.R. Larach 1973. An inversion in the atmosphere 
of Titan. Icarus 20, 437-446. 

Dierking, W. 1999. Quantitative roughness characterisation of geological surfaces and 
implications for radar signature analysis. Proceedings of the IEEE on Geoscience 
and Remote Sensing 37, n.5, 2397-2412. 

Elachi, C, E. Im, L.E. Roth, and C.L. Werner 1991. Cassini Titan radar mapper. 
Proceedings of the IEEE on Geoscience and Remote Sensing 79, n.6, 867-880. 

Eshleman, Von R., G.F. Lindal, and G.L. Tyler 1983. Is Titan wet or dry? Science 
221, 53-55. 

Frere, C, F. Raulin, G. Israel, and M. Cabane 1990. Microphysical modeling of Titan's 
aerosols: application to the in situ analysis. Adv. Space Res. 10, 159-163. 

37 



Fung, A.K., and H.J. Eom 1981. Multiple scattering and depolarization by a randomly 
rough Kirchhoff surface. Proceedings of the IEEE on Antennas Propagation AP-29, 
n.3, 463-471. 

Gibbard, S.E., B. Macintosh, D. Gavel, C.E. Max, I. de Pater, A.M. Ghez, E.F. Young, 
and CP. McKay 1999. Titan: high resolution speckle images from the Keck telescope. 
Icarus 139, 189-201. 

Griffith, C.A., T. Owen, G.A. Miller, and T.R. Geballe 1998. Transient clouds in 
Titan's lower stratosphere. Nature 395, 575-578. 

Griffith, C.A., J.L. Hall, and T.R. Geballe 2000. Detection of daily clouds on Titan. 
Science 290, 509-513. 

Guez, L., and P. Bruston, F. Raulin, and C. Regnaut 1997. Importance of phase 
changes in Titan's lower atmosphere. Tools for the study of nucleation. Planet. 
Space Sci. 45, n.6, 611-625. 

Hanel, B.A., and 15 colleagues 1981. Infrared observations of the Saturnian system 
from Voyager 1. Science 212, 192-200. 

Hansen, J.E., and L.D. Travis 1974. Light scattering in planetary atmospheres. Space 
Set. Rev. 16, 527-610. 

Ishimaru, A. 1997. Scattering and absorption of a wave by a single particle. In Wave 
Propagation and Scattering in Random Media (Donald G. Dudley, Series Editor) pp. 
9-40. Oxford University Press. 



38 



Khanna, R.K., M.A. Perera-Jarmer, and M.J. Ospina 1987. Vibrational infrared and 
Raman spectra of dicyanoacetylene. Spectrochim. Acta, 43A, 421-425. 

Khare, B.N., and C. Sagan 1973. Red clouds in reducing atmospheres. Icarus 20, 
311-321. 

Khare, B.N., C. Sagan, E.T. Arakawa, F. Suits, T.A. Callcott, and M.W. Williams 
1984. Optical constants of organic tholins produced in a simulated atmosphere: 
From soft X-ray to microwave frequencies. Icarus 60, 127-137. 

Lara, L.M., E. Lellouch, J.J. Lopez- Moreno, and R. Rodrigo 1996. Vertical distribution 
of Titan's atmospheric neutral constituents. J. Geophys. Res. 101, Issue E10, 23261- 
23284. 

Lara, L.M., E. Lellouch, J.J. Lopez-Moreno, and R. Rodrigo 1998. Erratum: Vertical 
distribution of Titan's atmospheric neutral constituents. J. Geophys. Res. 103, Issue 
Ell, 25775-25776. 

Lellouch, E., D.M. Hunten, G. Kockarts, and A. Coustenis 1990. Titan's thermospheric 
profile. Icarus 83, 308-324. 

Lorenz, R.D. 1993a. The life, death and afterlife of a raindrop on Titan. Planet. Space 
Set. 41, 647-655. 

Lorenz, R.D. 1993b. The surface of Titan in the context of ESA's Huygens probe. ESA 
J. 17, 275-292. 

Lorenz, R.D., and J.I. Lunine 1997. Titan's surface reviewed: the nature of bright and 
dark terrain. Planet. Space Sci. 45, 981-992. 

39 



Limine, J. I., D.J. Stevenson, and Y.L. Yung 1983. Ethane ocean on Titan. Science 
222, 1229-1230. 

McKay, CP. 1996. Elemental composition, solubility, and optical properties of Titan's 
organic haze. Planet. Space Sci. 44, 741-747. 

McKay, CP., J.B. Pollack, and R. Courtin 1989. The thermal structure of Titan's 
atmosphere. Icarus 80, 23-53. 

McKay, CP., A. Coustenis, R.E. Samuelson, M.T. Lemmon, R.D. Lorenz, M. Cabane, 
P. Rannou, and P. Drossart 2001. Physical properties of the organic aerosols and 
clouds on Titan. Planet. Space Sci. 49, 79-99. 

Muhleman, D.O., A.W. Grossman, B.J. Butler, and M.A. Slade 1990 . Radar reflec- 
tivity of Titan. Science 248, 975-980. 

Podolak, M., and R.E. Danielson 1977. Axel dust on Saturn and Titan. Icarus 30, 
479-492. 

Podolak, M., and L.P. Giver 1979. An inhomogeneous scattering model of Titan's 
atmosphere. Icarus 37, 361-376. 

Rages, K., and J.B. Pollack 1980. Titan aerosols: Optical properties and vertical 
distribution. Icarus 41, 119-130. 

Rages, K., and J.B. Pollack 1983. Vertical distribution of scattering hazes in Titan's 
upper atmosphere. Icarus 55, 50-62. 



40 



Rannou, P., M. Cabane, and E. Chassefiere 1993. Growth of aerosols in Titan's at- 
mosphere and related time scales: a stochastic approach. Geophys. Res. Lett. 20, 
967-970. 

Rannou, P., M. Cabane, E. Chassefiere, R. Botet, CP. McKay, and R. Courtin 1995. 
Titan's geometric albedo: role of the fractal structure of the aerosols. Icarus 118, 
355-372. 

Sagan, C, and W.R. Thompson 1984. Production and condensation of organic gases 
in the atmosphere of Titan. Icarus 59, 133-161. 

Samuelson, R.E., and L.A. Mayo 1997. Steady-state model for methane condensation 
in Titan's troposphere. Planet. Space Sci. 45, n.8, 949-958. 

Sen, A.D., V.G. Anicich, and T. Arakelian 1992. Dielectric constant of liquid alkanes 
and hydrocarbon mixtures. J. Phys. D 25, 516-521. 

Singh, S.P., and R.C. Miller 1979. Compressed-liquid dielectric constants and derived 
excess volumes for binary and ternary mixtures of N 2 , Ar, CH 4 and C 2 Hq. J. Chem. 
Thermodynamics 11, 395-409. 

Smith, B.A., and 25 colleagues 1981. Encounter with Saturn: Voyager 1 imaging 
science results. Science 212, 163-191. 

Smith, B.A., and 28 colleagues 1982. A new look at the Saturn system: The Voyager 
2 images. Science 215, 504-537. 

Smith, P.H., M.T. Lemmon, R.D. Lorenz, L.A. Sromovsky, J.J. Caldwell, and M.D. 
Allsion 1996. Titan's surface, revealed by HST imaging. Icarus 119, 336-349. 

41 



Thompson, W.D., and S.W. Squyres 1990. Titan and other icy satellites: Dielectric 
properties of constituent materials and implications for radar sounding. Icarus 86, 
336-354. 

Tomasko, M.G. 1980. Preliminary results of polarimetry and photometry of Titan at 
large phase angles from Pioneer 11. J. Geophys. Res. 85 A, 5937-5942. 

Tomasko, M.G., and P.H. Smith 1982. Photometry and polarimetry of Titan: Pioneer 
11 observations and their implications for aerosol properties. Icarus 51, 65-95. 

Toon, O.B., R.P. Turco, and J.B. Pollack 1980. A physical model of Titan's clouds. 
Icarus 43, 260-282. 

Toon, O.B., CP. McKay, R. Courtin, and T.P. Ackerman 1988. Methane rain on 
Titan. Icarus 75, 255-284. 

Toublanc, D. , J. P. Parisot, J. Brillet, D. Gautier, F. Raulin, and CP McKay 1995. 
Photochemical modeling of Titan's atmosphere. Icarus 113, 2-26. 

Ulaby, F.T., R.K. Moore, and A.K. Fung 1982. Microwave Remote Sensing: Active 
and Passive, Artech House, Noorwood MA. 

Ulamec, S., J. Badoz, and J. P. Lebreton 1992. Dielectric constant measurements in 
simulated Titan ocean liquids. In ESA SP-338 401-405. Proceedings Symposium on 
Titan. 

Veverka, J. 1973. Titan: Polarimetric evidence for an optically thick atmosphere. 
Icarus 18, 657-660. 



42 



Yung, Y.L., M. Allen, and J. P. Pinto 1984. Photochemistry of the atmosphere of Titan: 
Comparison between model and observations. ApJS 55, 465-506. 

Zellner, B. 1973. The polarization of Titan. Icarus 18, 661-664. 



43 



Table 1: Methane rain properties retrieved from modeling by Toon et al. (1988) and 
Courtin et al. (1995). 



Physical property Value 

Cloud extension 20 km 

Altitude of cloud top z = 30 — 40 km 

Mean raindrop radius between 50 \im and 3 mm (Toon only) 

Raindrop density between 1 and 3500 m~ 3 (Toon only) 

Optical thickness < 2 (far-infrared) 



Table 2: Set of samples supplied by simulations at LISA 



Synthesis conditions 



Volume of tholins deposited 



Set n°l 

(3 samples on glass tablets) 



V = 2.2 kV 
I = 80 mA 

gas mixture of N 2 /CH 4 (in proportion 98 : 2) 
p = 3 mbar 

deposit duration : 4h25 < 1.5 mm 3 

V = 2.6 kV 
I = 50 mA 

gas mixture oiN 2 /CH 4 (in proportion 98 : 2) 
p = 4 mbar 

deposit duration : 5h50 ~ 1.5 mm 3 

V = 2.6 kV 
I = 50 mA 

gas mixture of N2/CH4 (in proportion 98 : 2) 
p = 4 m6ar 

deposit duration : 17h00 ~ 8 mm 3 



Set n°2 V = 2.8 kV 

(1 glass tablet + 1 silica cylinder) 7 = 40 mA 

gas mixture of N2/CH4 (in proportion 98 : 2) 

p = 4 mfrar 

deposit duration : 6h45 < 8 mm 3 

~ 16 mm 3 



Set n°3 

(1 silica cylinder) 



deposit duration : 9h40 

V = 2.5 kV 
I = 30.5 mA 

gas mixture of N2/CH4 (in proportion 98 : 2) 
p = 1 m6ar 

deposit duration : 44h20 



2 mm 3 



Table 3: Tholin dielectric constants measured in resonant cavities. 



Cavity used Dielectric constant 

Set n°l 2.45 GHz e r = 2.26 - j0.12 

Set n°2 2.45 and 10 GHz e r = 2.2 - j0.002 

Set n°3 10 GHz e r = 2.03 - j0.022 



Table 4: Microwave dielectric constants of some major materials suspected to compose Titan's atmosphere and surface. 



Species 


Dielectric constant 


Frequency range of validity 


Source 




Pure liquid CH 4 (94 K) 


s r = 1.7- jl.5 x 10~ 2 


hundreds of MHz 


Thompson and Squyres 


(1990) 


Perfectly mixed 
CH 4 -C 2 H 6 liquid 
(around 100 K) 


s r = 1.8- j2 X 10~ 3 


GHz 


Singh (1979) 
Sen et al. (1992) 




Pure H 2 ice 
(between 90 and 273 K) 


£ ~ 3.1 

and very low loss tangent 


- 10 GHz 


Ulaby et al. (1982) 




Pure CH A -C 2 H 6 
ice (around 90 K) 


e - 2-2.4 

and very low loss tangent 


- 10 GHz 


Thompson and Squyres 


(1990) 



Table 5: Summary of the atmospheric models for Titan used in our simulations. 



Scenario label 



Description 



Physical properties 



Dry atmospheres 
A 

B 



only aerosols 

Homogeneous multi-layer model 

(with 0.1 /mi-radius particles in the main haze layer) 



Homogeneous multi-layer model 

(with 0.5 /xm-radius particles in the main haze layer) 



Tholin dielectric constant 
e = 2.2 ±0.1 
e" = 0.05 ±0.05 

idem 



C 
D 

Wet atmospheres 
E 



Inhomogeneous model taken from McKay et al. (1989) 
(see Fig. QKa)) 

Inhomogeneous model taken from Cabane et al. (1992) 
(see Fig. Q]-(a)) 

aerosols + condensation droplets 

Inhomogeneous model taken from Frere et al. (1990) 
including condensation at low alitude (below 100 km) 
(see Fig. EKa)) 

Homogeneous multi-layers model A 

+ methane rain layer between 30 and 10 km altitude 



Homogeneous multi-layers model B 

+ methane rain layer between 30 and 10 km altitude 

Inhomogeneous model C 

+ methane rain layer between 30 and 10 km altitude 
Inhomogeneous D 

+ methane rain layer between 30 and 10 km altitude 



idem 



idem 



Droplet dielectric constant 
(mixed hydrocarbons) 
e r = 1.8 -j2 x 10" 3 

Pure CHi rain layer 
(nominal case) 
extension: 10-30 km 
radius: 2 mm 

/ = 10~ 4 % 

e r = 1.7 -jl. 5 x 10~ 2 
idem 



idem 



idem 



FIGURE CAPTIONS 



Figure 1. (a) Inhomogeneous models of Titan's aerosol distribution inferred from mi- 
crophysical calculations : (1) McKay et al. (1989), (2) Frere et al. (1990), and (3) 
Cabane et al. (1992). (b) Schematic representation of differential condensation pro- 
cesses occurring in the 100 last kilometers of Titan's atmosphere as modeled by Frere 
et al. (1990). Differential condensation of hydrocarbons and nitriles control the size, 
density and chemical composition of the aerosol below 100 km. Aerosols are then of 
increasing size with decreasing altitude, reaching a maximum radius of ~ 900 \xm 
around 10 km altitude. Below this altitude, sublimation of liquids surrounding the 
aerosol leads to a reduction of the particle size until hypothetical complete sublimation 
occurs 2-3 km above the ground. 

Figure 2. Geometry and description of the multi-layer model we used for atmospheric 
transmission computation. 

Figure 3. (a): Attenuation versus particle radius and density for a 500 km thick 
homogeneous atmosphere with e r = 2.2 — j0.05; (b): 12.5 dB-threshold versus radius 
and density of dielectric spheres (e r = 2.2 — j0.05) for a 500 km thick homogeneous 
medium. The grayed region represents atmosphere models for which the radar is unable 
to see the surface of Titan. 

Figure 4. Attenuation caused by heterogeneous models (aerosol only and aerosol+condensation 



droplets). The description and properties of each scenario symbolized by letters in this 
figure can be found in Table 5. The minimum values for attenuation were calculated 
with e = 2.1 and e" = 0, and the maximum values with e = 2.1 and e" = 0.1. The 
mean attenuation was calculated with a value of e r = 2.2 — j0.05 contained within 
the two attenuation extrema symbolized by the errorbars. The comparison between 
"dry" models and "wet" ones emphasizes the drastic increase in attenuation due to 
rain occurrence. Errorbars are too small to be seen in the "wet" model case, due to the 
fact that aerosol attenuation - the only source of incertitude here- become negligible 
compared to the rain effect. The solid horizontal line shows the sensitivity threshold 
at -12.5 dB , and the dotted line represents the -4 dB threshold. 

Figure 5. (a): Attenuation for heterogeneous models AJj, BJj, CJj, and D$ with in- 
creasing droplet radius from no rain (0 mm) to the limit of 4.5 mm imposed by Lorenz 
(1993); (b): Attenuation for models AJj, BJJ, CJj, and D# with increasing condensa- 
tion efficiency on aerosols (factor /) varying between our nominal case (10~ 4 %), in 
agreement with Toon et al. (1988), and 0.1%. Errorbars on attenuation values due to 
uncertainties in the dielectric constant of tholin become insignificant beyond a certain 
value of r (2 mm) and are totally hidden by attenuation due to rain as soon as / 
is higher than 10~ 4 %. For (a) and (b), the dielectric constant of tholins was taken 
from the value we experimentally measured: e r = 2.2 ± 0.1 — j0.05 ± 0.05. The solid 
horizontal line shows the -12.5 dB sensitivity threshold and the dotted line the -4 dB 
threshold. 



Figure 6. Geometry of the single layer scattering problem. 

Figure 7. Compared radar backscattered power for a relatively smooth surface (a = 
0.5 cm and L = 25 cm, PO model) and a rough surface (a — 7 cm and L = 10 cm, GO 
model) covered by tholins (e r = 2.2 — j0.05), water ice (e r = 3.1 — jO.001) or liquid 
methane (e r = 1.7 — j0.015). The solid horizontal line shows the sensitivity threshold 
at -25 dB. 
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Figure 1: Rodriguez et al., Impact of Titan's aerosols on Cassini/Radar experiment 





Figure 3: Rodriguez et al., Impact of Titan's aerosols on Cassini/Radar experiment 
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Figure 4: Rodriguez et al., Impact of Titan's aerosols on Cassini/Radar experiment 
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Figure 5: Rodriguez et al., Impact of Titan's aerosols on Cassini/Radar experiment 
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Figure 7: Rodriguez et al., Impact of Titan's aerosols on Cassini/Radar experiment 



